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1— I ' Abstract 

o : 

^ ■ We calculate in this work the rates for the neutrino pair production by nucleon-nucleon 

bremsstrahlung taking into account the full contribution from a nuclear one-pion-exchange 
• 1—1 , 

^ ■ potential. It is shown that if the temperatures are low enough (T < 20MeV), the integra- 

H ; 

5r ! tion over the nuclear part can be done for the general case, ranging from the completely 
degenerate (D) to the non-degenerate (ND) regime. We find that the inclusion of the full 
nuclear contribution enhances the neutrino pair production by nn and pp bremsstrahlung 
by a factor of about two in both the D and ND limits when compared with previous 
calculations. This result may be relevant for the physical conditions of interest in the 
semitransparent regions near the neutrinosphere in type II supernovae, cooling of neutron 
stars and other astrophysical situations. 

1 Introduction 

Neutrino production processes in type II supernovae (SNII) and proto-neutron stars has 
received recently a considerable boost because of their importance in understanding of a 
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variety of interesting phenomena like the very explosion mechanism, aspects related to the 
SNII nucleosynthesis, the interpretation of the neutrino signal from SN, the thermalization 
and cooling mechanisms, among others. 

Several physical processes producing neutrinos have been investigated. The most im- 
portant ones are electron scattering, electron-positron annihilation, pion-nucleon and NN 
scattering, NN bremsstrahlung, plasmon- and photo-neutrino production. Their relative 
importance for different stages of the early proto-neutron star evolution (including the 
explosion itself) is a consequence of quite different functional dependences on p, T. An 
analogous statement can be carried over to the great deal of work done on the cooling of 
mature neutron stars 

One of the important questions in supernova physics is the actual post-bounce evolu- 
tion driven by neutrino emission, thought to be crucial for the long-term delayed shock 
revival driven by neutrinos (the Wilson mechanism). Neutrino production and thermal- 
ization were, for a long time, assumed to be dominated by pair annihilation and electron 
scattering respectively; at least for the i/^ and species (given that charged current re- 
actions are certainly the main source for t'e's). Recenty, however, increasingly accurate 
calculations revealed that the energy transfer rate in neutrino (z/^) - nucleon scattering 
may be up to one order of magnitude larger than the previous estimations, making this 
process competitive with the i/^-electron scattering as an equilibration mechanism @], 0|, 
0. Also, similar calculations for the total neutrino emissivity from NN bremsstrahlung 
indicate that this process may compete with e^e~ annihilation 0, [10|, [15] as a source 
for the and neutrinos. Bremsstrahlung rates have also been revisited for those non- 
degenerate to semidegenerate conditions, including the potentially important np channel. 
As a general feature of nucleon scattering, NN bremsstrahlung and URCA processes, it 
has proved quite difficult to treat the strong NN interaction matrix elements which are 
responsible for the v production together with the weak interactions. 

In the early calculations this interaction has been treated either by computing the over- 
lap integrals associated with the initial and final nucleon wave functions [|r^ or through 
the use of a Fermi liquid parametrization |T^. Later on, a NN potential based on the 
one-pion-exchange (OPE) approximation has been deduced by several authors 



0], |T3[, [0, p. However, in all these papers the full momentum transfer dependence 
of the OPE potential has been approximated (being, for instance replaced by the Fermi 
momentum) or even the whole square of the matrix element taken as a constant |0. 
However, for an accurate computation of these processes (and thus a realistic computation 
of the post-bounce evolution of type II supernovae) an accurate treatment of the nuclear 
potential is certainly required. 

We present in this work a method for treating explicitly the momentum dependence 
of the OPE potential involved in the calculation of the neutrino pair production from 
NN bremsstrahlung. Actually, the method is rather general and could be also used for 
the treatment of other processes involving NN potential. We show that in particular 
physical conditions, characterized by temperatures T < m\/m ~ 20MeV, the integral 
over the nuclear matrix element collapses, in a good approximation, to an integral which 
is independent of angles. The remaining part of the ME contribution, which depends only 
on the nucleon and neutrino energies, can be easily integrated numerically. Some results 
of this evaluation for non-degenerate and degenerate nucleons are presented. 

2 Calculations 

The total volumetric emissivity for the neutrino pair production by nucleon bremsstrahlung 
is given by 
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where 



F{f)=fMi -f3){i -h) (2) 

is the product of Fermi functions for the initial (1,2) and final (3,4) nucleons = 
( exp T +11 . In eq.(l) Pi, i = 1, 4 and q!/,qp are the nucleon and neutrino mo- 
menta, respectively; uj = uji+uj2 are the neutrino energies ; s is a symmetry factor taking 
into account the symmetry of identical particles (s = 1/4 for nn and pp, and s= 1 for the 
np channel) and Ei and /i, are the energies and chemical potentials of the nucleons. 



In the non-relativistic limit Ei ~ ''^ + ^ and one defines the chemical potential 
fi = fi—m. With the introduction the non-dimensional quantities [§[ y = f: ]Ui = 
we write fi = (exp"'^^' +1)"^. The degenerate (D) limit is achieved for y >> 1, while in 
the non-degenerate (ND) limit y << 1. 

In the OPE approximation the spin-summed NN matrix element (ME) for nn and pp 
interactions has been derived by several authors (see for instance 0) and reads 
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where qa = 1.26,/ ~ 1, A = QAG^Qa (i)^ k = pi - pg and 1 = pi - p4 are 
the nucleon direct and exchange transfer momenta, respectively. In the non-relativistic 
approximation one can also neglect the neutrino momenta in eq.(3), since it is always small 
compared to the nucleon ones. For the sake of completeness we quote the corresponding 
ME for the np interaction M 
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y k2 + m2 j 

It is convenient to perform the integrals in the center-of-mass coordinate system and 
define the new variables introduced by Brinkmann and Turner ^ 



Pi + P2 Pi - P2 

p+ = — 7s — i p = — 7s — ; 



P3c = P3 + P+; P4c = P4 + P+ 
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From the definition of tlie u variables above, one can easily deduce the following 
relations 



We now address the nn and pp expresion of the OPE potential and calculate the 
corresponding emissivity. After some algebra one can express the matrix element eq.(3) 
in terms of the scalar combinations = + P and = k^P as 

M - A _ m2 as(3ap + 2m^a,) + m^(7ap + 6m^a,) + 3ml \ ^ _ M'^rr^ /on 
~[ + 2<a.) + mlial + 2a,) + 2mla, + m| J " ^ 

Thus, the contribution of the nuclear ME can be sphtted into a (already calculated) 

constant term and a correction to be evaluated. It is worth mentioning that the main 

part of the nuclear potential contribution (the constant term) is just 3, i.e. the limit 

to which the expression eq.(3) converges when the pion mass is neglected compared to 

the nucleon momentum transfer. The a, and terms have the explicit expressions 

as — 2d{x + l/2z) ; a, = (P[x{x + z)sin^^ + 1/Az\ with d — 2mT, x — 2u3c and 

z = j;. Again, after some algebra the correction part is cast in the form 



- B,-B2{Cg + SgCOS(P)'' + Bs{Cg + SgCOS(l))^ ^' 

where Cg — cos^i ■ cos^c ', Sg — siwyi • siwyc and is the difference between the 
azithmutal angles corresponding to 71 and 7c (in spherical coordinates). 

The coefficients Ai and Bj are polynomials of degree eight in the product {rn^mT), 
depending only on variables x and z but otherwise independent of the angles. B^ is also 
binomial having a coefficient proportional to (2mT)^, while Bi and B2 contain terms 
proportional to m^, and m^. Therefore, whenever mT < (or equivalent, T < 
20 MeV), S3 can be neglected because its smallness compared to the other terms. This 
approximation can also be checked numerically. The results of integration over M^'^^ 
with and without the term Bs differ each other by a very small amount, (less than 2%) 
and greatly simplifies the calculation . The remaining integral over can be performed 
analytically yielding 
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The remaining integral over the other variables can be calculated by combining analytical 
and numerical techniques, for both the ND and D limits. 

3 Results 

We estimated first the total neutrino emissivity for the processes nnvv and ppvv since 
they involve the same contribution from the nuclear potential in both the ND and D 
cases. For comparison we normalize to the results of Thompson, Burrows and Horvath 
(2000) valid for a constant matrix element and arbitrary degeneracies. Denoting their 
result for the nn (or pp) neutrino emissivity with Q^" we find an nndenevaluation of their 
calculation of a factor of about two as compared with our result. More precisely we found 
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1.065 (11) 

in the ND limit. Furthermore, this difference is independent of y as long as the ND 
limit applies and T < m1/m . Care should be taken in this comparison since instead of 
sticking strictly to the high-momentum limit (a procedure that would have required the 
presence of the "3" prefactor), Thompson, Burrows and Horvath (2000) have introduced 
a fudge factor C, to embrace all those effects. In terms of that quantity we obtain ,^ ~ 2; 
or in other words, that the "true" bremsstrahlung emissivity is about 2/3 of its high- 
momentum limit. 

In the D limit, and adopting the same conditions of Flowers, Sutherland and Bond 
(1975) we obtain 



^ ^ = 1.0578 (12) 
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The neutrino emissivity for the case of the npuu process can be calculated analogously 
using an analogous procedure to obtain M™^^, but given that in general the chemical 
potentials of the neutrons and protons are different from each other this contribution is 



not so easily evaluated. We shall address the emissivity of this process and the general 
case for arbitrary degeneracy of all processes in a forthcoming paper. 



4 Conclusions 

In this work we developed a general method for calculating the total emissivity of the neu- 
trino pair production from NN bremsstrahlung, taken into account the full contribution of 
an OPE nuclear potential. We have showed that for particular physical conditions, char- 
acterized by T < 20 MeV, the multiple integral appearing in the emissivity formulae can 
be performed. Both in the ND and D regimes we found that the inclusion of the nuclear 
potential contribution produces neutrino emissivities which are ~ 2/3 of their respective 
high-momentum limits. The method for including the full contribution of a nuclear OPE 
potential in the calculation of the neutrino emissivites from NN bresstrahlung is rather 
general and allows the computation of all nn,pp and np processes in both ND and D 
limits and also it could be used to treat other proceeses of astrophysical interest, in which 
NN interaction is important. 

Some key issues still remain to be clarified for a full evaluation of the bremsstrahlung 
neutrino emissivity, the most important perhaps is the interplay between the correlations 
in the dense medium ([^, Q) and the momentum dependence of the nuclear potential. 
Other refinements (like the inclusion of multipion exchange) are not expected to be cru- 
cial, at least in the ND or mildly degenerate conditions corresponding to the semitranspar- 
ent post-shock conditions in a SNII, although they may be relevant in other situations. 
Another issue is the corresponence with the recently claimed reduction extracted in a 
model-independent analysis by Hanhart, Phillips and Reddy We have found that 



the inclusion of the full OPE potential indeed reduces the emission rates in both extreme 
limits, although a detailed comparison with that work has not been attempted here. 
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